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A piston is introduced into a spherical lune Casimir cavity turning it 
into two adjacent lunes separated by the (hemispherical) piston. On 
the basis of zeta-function regularisation, the vacuum energy of the ar- 
rangement is finite for conformal propagation in space-time. For even 
spheres this energy is independent of the angle of the lune. For odd 
dimensions it is shown that for all Neumann, or all Dirichlet, bound- 
ary conditions the piston is attracted or repelled by the nearest wall 
if d = 3, 7, . . . or if d = 1, 5, . . ., respectively. For hybrid N — D condi- 
tions these requirements are switched. If a mass is added, divergences 
arise which render the model suspect. The analysis, however, is rel- 
atively straightforward and involves the Barnes zeta function. The 
extension to finite temperatures is made and it is shown that for the 
3, 7, . . . series of odd spheres, the repulsion by the walls continues but 
that, above a certain temperature, the free energy acquires two minima 
symmetrically placed about the mid point. 
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1. Introduction. 



The Casimir piston, [1], has been the subject of considerable investigation, 
the object being to elucidate the nature of the Casimir effect after rendering some 
unwelcome infinities irrelevant from the outset. I refer only to the two recent works 
by Kirsten and Fulling, [2] and by Fucci and Kirsten, [3], where more extensive 
references and background information can be found. 

The basic situation reflects Casimir's original question concerning the global 
effect of introducing a partition into a cavity, Ai', so that Ai' becomes the union 
of two (joined) pieces, MU M*. 

Working with the free energy, F, this effect is expressed as the difference before 
and after constructing the partition: 



which was referred to as 'Casimir renormalisation' in [4]. As remarked in this 
reference, if the partition were flexible, or free, it would seek to adjust itself so as 
to minimise AF. 

The subtraction in (1) is not the same as that employed in [1], but has the 
same consequence. If the cavity, Ai', is considered to sit inside a larger enclosure 
(the universe), equation (1) for the internal quantities follows because the exterior 
contributions to the before and after free energies cancel, being independent of the 
partitioning. Furthermore, the last term in (1) is also independent of the partition 
and can be dropped when, say, evaluating the force on the partition. I thus regain 
the expression used in [1], and later works. 

At zero temperature, Lukosz, [5], uses a similar technique and I shall here 
mainly restrict myself to this case when F becomes E, the vacuum energy. Com- 
ments on finite temperature effects are made in a later section. 

There are a number of arrangements of this partitioning, that described in [1] 
being the most common, a transverse, movable flat plate (the piston) in a rectangu- 
lar box in flat space. By contrast, Fucci and Kirsten, [3], consider, as an example, 
a curved space situation where the manifold, Ai', is a (compressed) ball and the 
piston is an interior sphere enclosing a global monopole. 

Here, I wish to describe another example involving a curved space. I will not 
answer any outstanding questions, nor will the system be particularly realisable, 
but it has the advantage of a simple analysis, yielding closed formulae. One does 
not encounter the Epstein ^-functions arising in rectangular geometries, e.g. [4,6] 
and the many later discussions, which must be something of a novelty in this topic. 
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2. A curved piston 

My model takes space-time (the universe) as ultra static, of the form T x S 
where S is a <i-dimensional sphere whose azimuthal angle, 4>, has been unrolled so 
that < <p < oo. The cavity, Ai' , is obtained by choosing (d — 1) -hemispherical 
boundaries at, say, the azimuthal angles 4> = and 4> = a making it a <i-lune, S^, 
of angle a. I consider a scalar field in this cavity subject to Dirichlet or Neumann 
conditions at the boundary. The rest of the <i-sphere, another lune, of infinite angle, 
constitutes the 'outside' region. 2 These manifolds are singular. 

A piston 3 is now introduced in M' as another hemisphere located at the (vari- 
able) angle <j> = (3 < a, again with D or N conditions. This makes the manifolds 
M, M*, lunes of angles /3 and a — (3 respectively. There is nothing, in this model, 
to stop a becoming infinitely large, sending one wall to infinity. 

All one needs to implement the Casimir renormalisation, (1), is the vacuum 
energy on a D or N lune. 



3. Vacuum energy on lunes. Conformal coupling. 

It is a standard result, e.g. [4,6], that, if it is not infinite, the vacuum energy 
of a scalar field on an ultra static space time is given by, 

E=l( M (-l/2), (2) 
in terms of the (^-function on the spatial section, M., defined by 

where A n are the eigenvalues of the spatial part of the propagation operator. This 
part generally takes the form of the Laplacian, or improved Laplacian, on Ai. 

For simplicity, I assume that the scalar field is massless and is conformally 
coupled in d + 1 dimensions. This has the highly simplifying consequence of making 
the eigenvalues on the lune perfect squares, determined by, 

, 71" 

v /\ n = a + - ni +n 2 + ... + n d , (4) 



If one wished, the universe could be made finite by the introduction of a far boundary at some 
large value of or the unrolled sphere could be replaced by an ordinary one. All angles would 
then lie between and 2tv. 

The word 'piston' is, perhaps, slightly inaccurate. 
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where (3 is the angle of the lune, and the n i are non-negative integers. The constant 
a depends on the boundary conditions, Neumann conditions corresponding to a = 
a N = (d — l)/2 and Dirichlet to a = a D = /3/n + (d — l)/2.. 

These eigenvalues follow by standard separation of variables. Gromes, [7] , gives 
the analysis for the two-lune (see also, earlier, Pockels, [8] pp. 100-109). The higher 
lune result reflects the nesting of their metrics. The resulting ^-function from (3) 
is a Barnes (^-function. 

The situation is exactly the same as that encountered in [9] except that there 
the lune is a special one with (3 = Tc/q, q G Z. The vacuum energies can then 
be read off from the expressions in [9] because the particular analysis extends to 
arbitrary f3. 

For completeness, I give some appropriate formulae. The (^-functions (for con- 
formal coupling in d + 1 dimensions) are 

C JV (a) = Cd(2a,(d-l)/2|7r/ j 9,l) 

Cd(s) = Cd(2s, tt//3 +{d- l)/2 | 7T//3, 1) 
in terms of the Barnes (function, £ d . Looking at (2), one needs the value of the 
Barnes function at — 1, and Barnes himself, [10], gives this as a generalised Bernoulli 
polynomial. The expression for the vacuum energies on lunes is, 

E{fi) = -(T) d+ \^ +iy B% 1 ((d-l)/2 | n/f},l). (6) 

The upper sign is for Neumann and the lower one for Dirichlet conditions. Inciden- 
tally, adding the iV and D values gives the energy on a periodic lune of angle 2(3 
and it can be seen that this vanishes for even d, as is well known for the full sphere 
(for conformal coupling). 

Explicit computation of the polynomial in (6) reveals that E(/3) is independent 
of (3 for even d. This can be shown analytically as follows. 

By applying Barnes' derivative formula, 



I find 



~Bl d \a\ q,u) = -1 B[ d+1 \a + q\ q, q, W ) , 
oqq q z 

= -(-l) fc 1 Bl d+1) (d - 1 - a + q\ q, q, u>) , 



§~ q l q BflM ~ 1)/2| q, 1) = -1 B%?\(d - l)/2 + q\ q, q, 1) , 

= (-l) d ^B^\ 1 \(d-l)/2 + q\q,q,l) 



(7) 



(8) 
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which therefore vanishes if d is even. 

This result means of course that there is no force (torque) on the internal 
partition for even d. 

Turning now to the more interesting case of odd d, I give two examples of the 
rational functions for the internal vacuum energy on an odd D or N lune of angle 

P = */q, 

2q 6 + Uq 4 + 77q 2 - 62 , r 

= , d = 5 . 

120960g 

The Casimir renormalised internal energy of the partitioned manifold is 

E((3, a) = E(fi) + E(a — 0) — E(a) . (9) 

As a function of (3, E(j3,a) is symmetric about the middle point, (3 = a/2, where, 
therefore, it has an extremum. If a minimum (stability) , the partition is repelled 
by the nearest wall, while if a maximum (instability), it is attracted. Repulsion is 
more interesting. 

Calculation shows that (9) has a minimum for d = 3, 7, 11, . . . and a maximum 
for d = 1,5, 9, 13, .... I obtain this behaviour analytically in Appendix 1. A roughly 
similar alternating behaviour is found by Fucci and Kirsten, [3] , in a quite different 
setup, for Dirichlet conditions (but not, curiously, for Neumann ones). 

Sending one wall to infinity a — > oo while keeping (3 finite does not alter these 
conclusions since one is working on the left-hand side of the energy (potential) 
curve. 



4. Hybrid boundary conditions 

One can ring the changes on boundary conditions. For example, still retaining 
symmetry, Neumann conditions could be enforced at the outer walls, with Dirichlet 
at the piston. 

The energy of an ND lune is given in terms of that on a DD, or iViV, one, (6), 
by (e.g. [11]), 

E ND (/3) = E(2/3) - E(/3) , 

and, following the analysis through, the dimension dependence of a maximum or a 
minimum at the mid point is switched so that, this time, one has a maximum for 
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d = 3, 7, 11, . . . and a minimum for d = 1,5, 9, 13, . . .. For even d, the ND energy is 
zero. 

In dimension d = 1, the consequent repulsion agrees with an old result, dis- 
cussed and extended by Fulling, Kaplan and Wilson, [12]. 



5. Finite temperature 

Another consequence of the simple form of the eigenvalues is that the finite 
temperature quantities are given as simple series, easily computed. For example 
the free energy at temperature T is, [9], on the lune of angle (3, 

F(T,P) = E{fi) -( coth ^ ± l^osech^ . (10) 

m=l 

As T increases, any minimum at the mid point will, at some point, become a 
maximum. This is best seen from the high temperature expansion of the free energy 
given in [4] and [13]. When computing the subtracted quantity, (1), the first two 
terms in this series, which are proportional to the measure of the domain concerned 
and of its boundary, respectively, cancel or are constant. The third term is the 
active one. It is proportional to the heat-kernel coefficient, C x , and equals, 
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The subtracted quantity, (1), constructed from (11) always has a maximum at 
(3 = 7r, whatever the dimension, d. (I am again setting a = 2n for convenience.) 
This is a limiting behaviour, and to obtain the precise dependency a numerical 
evaluation is required. This would yield the value of the temperature at which any 
central minimum just turns into a maximum 4 . The finite temperature correction 
(the second term in (10)) is always negative and has a maximum at the mid point. 
The value of F is a competition between the two parts of (10). 

As the walls are approached, the zero temperature part of the free energy, (10), 
dominates, again giving a repulsion for the odd dimensions, 3, 7, . . .. The graph of 



4 It is really the value of the dimensionless quantity, aT, that is determined, where a is the radius 
of the lune. For a = 2n the value of this transition temperature is approximately 0.455 for 
Neumann conditions in three dimensions. 
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the free energy against the angle of the piston is that of a double well potential 
so that the piston, this time, would preferentially settle in one of the two minima 
symmetric about the mid point, a somewhat curious situation, if it could be realised. 
For T = 0.6 these minima are very approximately at (3 = tc/2 and (3 = 3n/2. In 
even dimensions, a finite temperature provides an attraction to the walls. 

For the record, I give here the first three terms in the high temperature expan- 
sion of the free energy on a 3-lune of angle a, 

F(T, a) ~ - ™ - Cb(3) T 3 _ t!l (l _ « W , 
1 ' ' 90 4 36 \a n) 

which is useful when checking the exact form, (10). 
6. Conclusion 

In the 'physical' dimension, d = 3, the hemispherical piston arrangement is 
stable and would oscillate about the midpoint, if the temperature were below the 
transition value. For higher temperatures, the equilibrium position would shift to 
one of two minima symmetrically placed with respect to the midpoint. 

By contrast, for d = 5 say, the system is unstable, and the piston would collapse 
to one of the walls, whatever the temperature. 

It seems difficult to maintain that the system analysed in this paper is a prac- 
tical one. Rather, it should be considered as a mathematical fiction whose value 
lies in a relatively easy analysis. 

It has to be admitted that this simplicity is, to some extent, accidental, being a 
result of the conformal propagation equations which allow ( M (—l/2) to be finite. If 
these are modified, say by adding a mass, then divergences associated with conical 
singularities appear which are not removed by the subtraction, (9), nor, worse, are 
not independent of the position of the piston. (See Appendix 2.) That this does 
not occur in the rectangular piston system is a consequence of the fact that there 
the constituents are always perpendicular. If the piston were allowed to slant, even 
slightly, not only would the force change but it would acquire a divergence, which 
raises questions of principle. 

The system described in the present paper is a curved space analogue of the 
traditional wedge geometry investigated over the years. Most discussions involve the 
energy density, but there are some of the total Casimir energy in a bounded wedge. 
Recent calculations that make relevant remarks on divergences in the light of various 
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boundary conditions are Milton, Wagner and Kirsten, [14], Brevik, Ellingsen and 
Milton, [15] and Ellingsen, Brevik and Milton, [16]. 

In order to more closely mimic the wedge bounded by a circular arc in flat 
space, it is possible to impose boundary conditions at the equator of the <i-sphere 
which produces half a <i-lune. The eigenvalues are a simple modification of (4) and 
details will be exposed at another time. 



Appendix 1. Alternating extrema 

I prove that the vacuum energy, (9), is a maximum or a minimum at the mid 
point depending on the sign (— l)( d + 1 )/ 2 , for odd d. I do this by computing the 
second derivative of E(f3, a) with respect to j3 at f3 = a/2. 

There is nothing special about any particular value of the outer wall angle, a, 
which I now set equal to 2ir. The mid point is then at j3 = it, i.e. q = 1, its most 
convenient value. 

At the mid point, the first derivatives of the first two terms in (9) cancel, while 
the second derivatives are equal, and add. Hence it is sufficient to evaluate, using 

(6) 

(**&((d-l)/2|i 1)) (12) 

at x = 1 where x = (3/tt. In order to more directly employ Barnes' formula, (7), I 
use the variable q = 1/x and therefore seek, 

^ ^q \ B » d){a 1 ^ W) = q2 §-q Bijl+1) (a + Q 1 ^ ^ ^ ' (13) 

for the special values, v = d + 1, a = (d — l)/2, q = 1 and uj = 1. 

The three dependencies on q on the right-hand side of (13) are treated sepa- 
rately, and, using the symmetry of the parameters (q, q' , u>) in q and q', I find the 
general formula, 

^Bi d+1 \a + q\q,q,u>) = 



dx 



2 



+ ^-Bl d + 1 \a\q 1 q^) 

a — a + q 0(1 
q'=q 

-B^ (a + 2q\ q, q, q, w) + - (a + q\q,q,u>) 
q q 



a=a+q (14) 



+ uB { J L X 1 \a + q I q,q,u) 
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where I have used (7) again, and a standard identity for Bernoulli polynomials. 

Inserting the special values, it is seen, firstly, that the Bernoulli polynomials 
have all their parameters equal to unity, the most extensively discussed case. Sec- 
ondly, for d odd, the arguments a + 2q = (d + 3)/2 and a + q = (d+l)/2 are positive 
integers and, furthermore, are such as to make the first and last terms on the final 
line vanish. This is because of the product structure (e.g. [17] p. 186, [18], §8), 

B { n n+1 \x\ l) = (x- l)(x - 2) . . . (x - n) , 

which always possesses a vanishing factor. This leaves the second term, whose 
contribution to the second derivative at the mid point is 

4< + 1 1) ((^+l)/2|l) 
and it is the sign of this quantity that I require. I use the representation 

rx+l 

B { n n \x\ 1) = / dt(t- l)(t - 2) ... (t - n) 

J X 

and obtain, after translating the integration variable, 

It is easy to see that in the integration range, —1/2 < t < 1/2, the integrand is 
positive for d = 3,7, 11, . . . and negative for d = 5, 9, 13, . . .. This means that the 
vacuum energy is a minimum at the mid point in the former cases and a maximum 
in the latter, which was to be proved. 

Appendix 2. Divergences for non— conformal coupling 

To illustrate the emergence and character of the divergences, it is sufficient to 
add a small mass term, dm 2 , to the conformal eigenvalues A n , (4), and employ a 
binomial, perturbation expansion to give the resulting change in the (^-function, 

5( d (s, a | w) = -5m 2 s Q(2s + 2, a | u>) (15) 

which, with the appropriate a and u>, should be added to the conformal ^-functions, 
(5). The correction, (15), possesses a pole at the evaluation point s = —1/2 with 
residue proportional to the Bernoulli polynomial, 

Bf} 1 ((d-l)/2\n/P,l), 
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examples of which are the well known forms, associated with conical singularities, 



(7r 2 -/3 2 )(7r 2 + ll/? 2 ) 

30^ ' d " 5 ' 

and are, up to a factor, the heat-kernel coefficients, C( d _ 1 - )/ / 2 ? on t ne hme. For even 
d, the Bernoulli polynomial is always proportional to just n/(3. 

These results show that the /3-dependence of the divergence does not disappear 
from the vacuum energy combination (9). This presents us with an old problem of 
interpretation which is not resolved here. For this reason I do not take (15) any 
further. (Even though the number of pole terms is finite for any given d, when the 
higher perturbation terms are evaluated.) 
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